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(57)Abstract: 

PURPOSE: To efficiently prepare new mesh data in a 
short time at the time of changing the compression 
states of mesh by changing a part of component size to 
be the reference of compression control after preparing 
the mesh data and changing the mesh data by smoothing 
to the size after the change. 

CONSTITUTION: A command analysis part 101 analyzes 
commands from a user and starts a mesh preparation 
part 102, a compression correction part 103 and a l 
display part 104. The mesh preparation part 102 refers 
to shape data stored in a shape data storage part 105 
corresponding to a mesh preparation control command 
from the command analysis part 101, divides the shape 
to the mesh, prepares the mesh data and stores them in 
a mesh data storage part 106. The compression 
correction part 1 03 refers to the shape data 
corresponding to a compression correction control 
command from the command analysis part 101 and 
updates the mesh data of the mesh data storage part 

106 based on a part of components to be the reference of compression state change included 
in the command and the specified size after the change of the component. 
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Partial English Translation of Japanese Laid-Open Patent 
Application NO.H6-274573 

[0018] FIG. 1 is a diagram showing a configuration according 
to an embodiment of the present invention. The configuration 
comprises a command analyzing unit 101, a mesh forming unit 102, 
a roughness/denseness correction unit 103, a display unit 104, 
a shape data storing unit 105, and a mesh data storing unit 106. 
[0019] In this configuration, the command analyzing unit 101 
analyzes an instruction from a user to operate the mesh forming 
unit 102, the roughness/denseness correction unit 103, and the 
display unit 104. 

[0020] The mesh forming unit 102 refers to shape data stored 
in the shape data storing unit 105 according to a mesh formation 
control instruction from the command analyzing unit 101 to 
divide the shape into a mesh, thereby forming mesh data. The 
formed mesh data is stored in the mesh data storing unit 106. 
[0021] The roughness/denseness correction unit 103 refers to 
shape data according to a roughness/denseness correction 
control instruction from the command analyzing unit 101 to 
update the mesh data stored in the mesh data storing unit 106 
on the basis of an element serving as a reference of a change 
in roughness/denseness state included in the 
roughness/denseness correction control instruction and a 
designated size of the element obtained after the 
roughness/denseness state is changed. 

[0022] The display unit 104 converts the mesh data into a 
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graphic according to a display control instruction from the 
command analyzing unit 101 and display the graphic. 
[0023] The details of the shape data and the mesh data will 
be described below. 

[0024] The shape data is constituted by an apex table 201 in 
FIG. 2, a side table 301 in FIG. 3, and a surface table 401 in 
FIG. 4. 

[0025] The contents of the tables 201 , 301 , and 401 will be 
described below by using a three-dimensional cubic bl as shown 
in FIG. 5 as an example. 

[0026] In the apex table 201 , for respective identification 
codes (apex IDs) vl to vlO, coordinate values (xl, yl, zl) to 
(xlO, ylO, zlO) of the apexes vl to vlO are stored. 
[00.27] In the side table 301, for respective identification 
codes (side IDs) 11 to 115, apex IDs vl to vlO of both the ends 
(start points and end points) of sides 11 to 115 and curve data 
ell to cll5 are stored. In this case, the curve data indicates 
the radius and central coordinates of an arc, the control point 
coordinate string of a free curve, and the like. 
[ 0028 ] In the surface table 401 , as shown in an exploded diagram 
in FIG. 5(b) , for respective identification codes (surface IDs) 
si to s7 of surfaces, the number of sides constituting the 
surfaces si to s7, side IDs, and curved surface data csl to cs7. 
The curved surface data is expressed by coordinates or the like 
of a control point of a curved surface patch. 
[0029] When a user instructs the command analyzing unit 101 
to divide the three-dimensional cubic bl shown in FIG. 5(a) into 
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a mesh, the mesh forming unit 102 operates to generate a mesh 
which is divided as shown in FIG. 6. 

[0030] The mesh in FIG. 6 is expressed by data in a node table 
7 01 in FIG. 7, an element table 802 in FIG. 8, a node coupling 
table 901 in FIG. 9, and a shape-element-node-corresponding 
table 1011 in FIG. 10. 

[0031] The node table 701 of these tables is constituted by 
coordinate values (xl, yl, zl) to (x80 f y80, z80) classified 
by node numbers nl to n8 0 of the mesh. In this case, the node 
numbers nl to n80 are sequentially added to apexes, sides, 
surfaces, and the inside of a cubic, and the node numbers nl 
to n8 0 are stored in this order. On a side, node numbers are 
sequentially stored from the start point of the side to the end 
point of the side. 

[0032] An element table 801 is constituted by node numbers 
constituting elements indicated by element numbers el to e3 6 
of a hexahedron partitioned by a mesh. 

[0033] For example, the element number el is constituted by 
node numbers nl, nil, n40, nl9, n29, n51, n70, and n65, 
[0034 ] The node coupling table 901 is constituted by the number 
of nodes at which the elements of the node numbers nl to n80 
are connected by element sides and the node numbers of the nodes . 
[0035] In addition, the shape-element-node-corresponding 
table 1011 is constituted by the number of nodes generated on 
shape elements ( apexes , sides , surfaces, and cubics) and 
addresses (subscripts of an alignment) in the node table 701 
in which corresponding node number strings are stored. In this 
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case, a shape element ID corresponds to each ID (vl to vlO, 11 
to 115 , si to s7, bl) of the shape data shown in FIGS. 2 to 4 . 
An arrangement of apexes in the table is equal to an arrangement 
of corresponding nodes in the node table 701. The same applies 
to the sides and the surfaces. 

[0036] When a user instructs the command analyzing unit 101 
to multiply the size of an element connected to, e.g. , the node 
nl6 by a for the mesh in FIG. 6, the roughness/denseness 
correction unit 103 and the display unit 104 operates. 
[0037] For example, a is changed at intervals of 0.1 in the 
range of 0.7 to 0.3, and an instruction is made to concentrate 
elements on the node nl6 . In this case, the roughness/denseness 
correction unit 103 executes roughness/denseness correction 
at respective a values, and the display unit 104 sequentially 
displays results. 

[0038] The process performed by the roughness/denseness 
correction unit 103 will be described below with reference to 
FIG. 11. 

[ 0039 ] An element the size of which is multiplied by a and nodes 
constituting the element are extracted (step 1101 ) . In the case 
in FIG. 6, elements el 0 and ell having the node nl 6 as constituent 
nodes are retrieved from the element table 801. At this time, 
a node string N = {nl5, nl6, nl7, n39, n41, n43, n59, n60, n61, 
n69, n71, n73 Constituting the elements elO and ell are 
simultaneously obtained. The respective nodes of the node 
string N are nodes for controlling roughness/denseness 
correction. 
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[0040] Shape elements to which the respective nodes of the node 
string N are extracted (step 1102). For this purpose, node 
table addresses of the shape-element-node-corresponding table 
1011 are used, and node strings on all the sides are retrieved 
from the node table 7 01 to check whether the nodes of the node 
string N exist on the node strings on the sides or not. 
[0041] Furthermore, the nodes of all the surfaces will be 
checked. In the examples in FIGS. 5 and 6, of the node strings 
N, a node string NL1 = {nl5, nl6, nl7} belong to a side 14, a 
node string NS1 = <n3 9, n41, n43} belongs to a surface si, and 
a node string NS2 = {n59, n60, and n61} belongs to a surface 
s6. 

[0042] Node coordinates on the extracted sides are updated 
(step 1103). In this updating operation, with respect to the 
sides extracted in step 1102, a common ratio is calculated from 
the lengths of the sides of the element multiplied by a and 
the number of divisions, and the sides are subdivided such that 
the lengths of the element sides constitute a geometric 
progression, so that the node coordinates of the node table 701 
are updated. 

[0043] Since the size of an element connected to the node nl6 
is multiplied by a, with respect to a side 14, a length nl5 
- nl6 and a length nl6 - nl7 of the element sides are multiplied 
by a to obtain ell and el2. A length from the node nl6 to a 
start point v4 of the side 14 and a length from the node nl6 
to an end point v5 of the side 14 are calculated by using curve 
data of a side table 402. The lengths are represented by til 
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and tl2, respectively. 

[0044] Since the number of divisions from the node nl6 to the 
apex v4 is "2", a common ratio r obtained when the node nl6 and 
the apex v4 are used as a start point and an end point, 
respectively, can be calculated by: 
(r 2 - l)/(r - 1) = tll/ell 

For this reason, this equation is calculated to obtain the 
following result: 
r = tll/ell - 1. 

[0045] By using the common ratio r calculated as described 
above and the curve data in the side table 301, the coordinates 
of the node nl5 will be calculated- Similarly, a common ratio 
r obtained when the node nl6 and an apex v5 are set as a start 
point and an end point, respectively, and the coordinates of 
a node nl7 is calculated by using the curve data in the side 
table 301. 

[0046] Node coordinates node coordinates on the respective 
extracted surfaces are updated (step 1104 ) . This process means 
that node coordinates in the node table 7 01 with respect to the 
surfaces extracted in step 1102 are updated. The details of 
the process is shown in FIG. 12. 

[0047] A local coordinate system is set first (step 1201). 
This process means that a local coordinate plane for smoothing 
node coordinates is calculated. A start point of the first side 
of the surface is set as an origin o, and a vector extending 
from the start point of the first side to the end point of the 
first side is set as an x axis. An average normal vector of 
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the surface is set as a z axis, and an orthogonal coordinate 
system o-x-y-z is set such that the vector product of the z axis 
and the x axis is equal to a y axis, 

[0048] A coordinate transformation matrix is calculated (step 
1202)- In this case, a transformation matrix GL used in 
transformation from the whole coordinate system into a local 
coordinate system and an inverse matrix LG thereof are 
calculated. 

[0049] Coordinate transformation from the whole coordinate 
system into the local coordinate system is performed (step 12 03 ) . 
Data subjected to coordinate transformation includes the 
coordinate values of nodes on a surface, nodes on apexes and 
sides constituting the surface, control point of curved surface 
data. 

[0050] The apexes and the sides constituting the surface can 
be retrieved from the surface table 401 and the side table 301, 
the coordinates of the nodes on the surface, the sides, and the 
apexes can be extracted from the node table 1011 with reference 
to node table addresses of a shape-element-node-corresponding 
table 504. 

[ 0051 ] The control point coordinates of the curved surface data 
are obtained from the surface table 401. At this time, local 
z-coordinate values of all the transformed nodes are set at "0" . 
[0052] Of the constituent nodes of an element the size of which 
is multiplied by a, the local coordinates of the nodes on the 
corresponding surface are calculated ( step 1204 ) . At this time, 
the local coordinates of the nodes are changed such that the 
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distance between the nodes extracted in step 1102 and the center 
of the element the size of which is changed is a times. 
[0053] For example, when a surface s6 is used, the nodes of 
the node string NS2 may be moved such that distances between 
the nodes of the node string NS2 and the node nl6 on the surf ace 
of a local coordinate system are a times . 

[0054] The local coordinates of the nodes on the surface are 
smoothed (step 1205) . 

[ 0055 ] The process of smoothing the node coordinates means that 
the coordinate values of the respective nodes are repeatedly 
updated by an average value of the coordinate values of a 
plurality of nodes connected to the corresponding nodes to 
smooth the node coordinates. In this case, the connected nodes 
are registered in the node coupling table 901. 
[0056] Only nodes on the surface except for the nodes the 
coordinates of which are calculated in step 1204 are smoothed. 
The nodes on the apexes and the sides the coordinate values of 
which have been changed are not smoothed. For this reason, upon 
completion of this process, on the x-y plane of the local 
coordinate system, the nodes on the surface are arranged such 
that the element the size of which is multiplied by a can be 
smoothly attracted. 

[0057] When the corresponding surface is curved, nodes on the 
surface are projected on the curved surface (step 1206). The 
projected points are defined as intersection points between a 
straight line perpendicular to the x-y plane of the local 
coordinate system passing through the nodes on the surface and 
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curved surface patch expressed by the curved surface data 
subjected to coordinate transformation in step 1203. This 
process is performed to all the nodes of the surface to calculate 
z-coordinates in the local coordinate system. 
[0058] Coordinate transformation from the local coordinate 
system into the whole coordinate system is performed to all the 
nodes on the surface (step 1207). 

[0059] The above ends the updating of the coordinates of the 
nodes on the surface (step 1104). 

[0060] Finally , node coordinates in a cubic bl are updated 
(step 1105) . This process means that internal node coordinates 
of the cubic bl are extracted from the shape-element-node- 
corresponding table 1011 and the node table 7 01 to perform 
smoothing. 

[0061] The above is the process of the roughness/denseness 
correction unit 103. According to the process, 

roughness/denseness correction at certain a, and mesh data is 
updated. 

[0062] Upon completion of updating of mesh data, the display 
unit 104 operates in response to an instruction from the command 
analyzing unit 101, the display unit 104 displays updated mesh 
data . 

[0063] A result of a change in roughness/denseness state when 
a = 0.5 is shown in FIG. 13. 

[0064] After the result of the change is displayed by the 
display unit 104, the command analyzing unit 101 asks a user 
whether roughness/denseness correction is ended or not. When 
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a user wants to end the roughness/denseness correction, the user 
may input an end instruction in the command analyzing unit 101. 
When the user wants to change a to further execute the 
roughness/denseness correction/ the user may input a continuous 
instruction. 

[0056] Therefore , when a user wants to change a 
roughness/denseness state, the user designates an element 
serving as a reference of a change in roughness/denseness state 
and the changed size of the element. Only the related node 
coordinates on the sides and the surfaces and the node 
coordinates in the cubic are updated in the node table 701. 
[0066] Since division of a mesh is performed by a smoothing 
process, data obtained after the roughness/denseness state is 
changed can be obtained at a speed higher than that in a case 
in which a mesh is divided from the first by changing division 
ratios of sides. 

[0067] In addition, since an element serving as a reference 
of a change in roughness/denseness state can be designated on 
not only sides but also surfaces or cubics, the 
roughness/denseness state can be easily changed in the surfaces 
or the cubics . 

[0068] When the command analyzing unit 101 transmits a command 
to the roughness/denseness correction unit 103 such that the 
size of an element serving as a reference of a change in 
roughness/denseness state, the roughness/denseness states can 
be interactively changed while checking the 
roughness/denseness states at the display unit 104 in 
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respective stages . 

[0069] In this case, the gradual change in size may be manually 
performed or automatically performed such that the size is 
changed at predetermined change intervals. 

[0070] The present invention can be applied to not only the 
three-dimensional shape but also a two-dimensional shape. 
When the present invention is applied to the two-dimensional 
shape, the process in step 1105 in FIG. 11 is not necessary. 
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